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We report measurements of the beam spin asymmetry in Deeply Virtual Compton Scattering 
(DVCS) at an electron beam energy of 4.8 GeV using the CLAS detector at the Thomas Jefferson 
National Accelerator Facility. The DVCS beam spin asymmetry has been measured in a wide range 
of kinematics, 1.0 (GeV/c)^ < < 2.8 (GeV/c)^ 0.12 < xb < 0.48, and 0.1 {GeV/cf < -t < 0.8 
(GeV/c)^, using the reaction ep epX. The number of H(e,e'7p) and H(e,e'7r°p) events are 
separated in each {Q^,XB,t) bin by a fit to the line shape of the H(e,e'p)X distribution. The 
validity of the method was studied in detail using experimental and simulated data. It was shown 
that with the achieved missing mass squared resolution and the available statistics, the separation 
of DVCS-BH and n'^ events can reliably be done with less than 5% uncertainty. The Q^- and t- 
dependences of the sin moments of the asymmetry are extracted and compared with theoretical 
calculations. 

PACS numbers: 



PACS numbers: 13.60.Z, 14.20.Dh, 24.85. +p 



I. INTRODUCTION 
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Hard scattering processes play an important role in 
the understanding of the quark and gluon structure of 
hadrons. The important feature of hard reactions is 
the possibility of separating the perturbative (short dis- 
tance) and non-perturbative (long distance) parts of the 
interaction. This so-called factorization property has 
been successfully used in inclusive measurements {e.g. 
in Deep Inelastic Scattering (DIS) of leptons) to study 
the internal structure of the nucleon. Until recently, very 
few exclusive processes could be treated in the framework 
of perturbative QCD (pQCD) and compared to experi- 
mental data (typical examples are the 71^77* transition 
form factor T| and the elastic form factors of the pion 
[2] and the nucleon [3||). The recently developed formal- 
ism of a QCD description of Deeply Virtual Compton 
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Scattering (DVCS) 0, i] and Deeply Exclusive Meson 
Production provides a framework in which the am- 
plitudes of these processes can be factorized into a hard- 
scattering part (exactly calculable in pQCD) and a non- 
perturbative nucleon structure part that can be param- 
eterized at the amplitude level by means of Generalized 
Parton Distributions (GPDs). The GPDs contain in- 
formation on quark/ antiquark correlations, particularly 
the correlation of their transverse spatial and longitudi- 
nal momentum distributions, and on the quark angular 
momentum Q. They provide a unifying picture for an 
entire set of fundamental quantities containing informa- 
tion on the hadronic structure, such as nucleon form fac- 
tors (which are related to matrix elements of vector and 
axial vector currents), polarized and unpolarized parton 
distributions, and contributions to the spin of the nu- 
cleon due to orbital excitations of quarks and gluons. 

There are four chiral-even GPDs, denoted H'', H'^, 
i?', and i?', which depend on the kinematical variables 
X, 5, and t. They correspond to the amplitude for remov- 
ing a quark with momentum fraction x + ^ and restor- 
ing it with momentum fraction x — ^. The light-cone 
momentum fraction x is defined by fc+ ~ xP~^ , where 
k is the quark loop momentum and P is the average 
nucleon momentum {P — {p' +p)/2, where p and p' 
are the initial and final state nucleon four-momenta, 
respectively). ^ is the generalized Bjorken variable, 
( = QV(4g ■ P) 2:5/(2 - xb) as t/Q^ 0, where 
q = k — k' and — —q^- k and k' are the initial and 
final electron momenta. 

The Mandelstam variable t — = {p' — pY is the 
Lorentz-invariant four momentum transfer squared to 
the target. In the forward limit, t ^ 0, the GPDs H 
and H reduce to the quark density distributions q{x) and 
quark helicity distributions Aq(a;) obtained from DIS: 



ular quark flavor (for any f): 

.+1 

j dxHi{x,^,t) ^ F^{t) 



dxE^x,^,t) = F^{t) 



dxH^{x,C,t)^g\it) 



+1 



dxE'i{x,i,t) = h\{t), 



(3) 



(4) 



(5) 



(6) 



where FKt) and -F||(i) represent the elastic Dirac and 
Pauli form factors, respectively, for the quark flavor q 



in the nucleon. 
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is the axial-vector form factor of the 



nucleon, and h\ is the pseudo-scalar form factor. 

Deeply Virtual Compton Scattering (DVCS), ep — > 
ep7, is the simplest reaction to access GPDs experimen- 
tally. DVCS has the same final state as the Bethe-Heitler 
(BH) process, where a photon is emitted by the incom- 
ing or outgoing electron (Fig. [T]). The differential cross 
section for the process ep — > ep7 can be written as Q : 



d<7 



a'^XBU 



dxBdydtd<j) SttQ'^^/T 



\T\ 



(7) 



Here the amplitude of the production of a real photon, 
r, is the sum of the DVCS {Tdvcs) and BH (Tbh) 
amplitudes, given as: 



BH\ 



IT, 



DVCS\ 



(8) 



F9(x,0,0) = 
F«(x,0,0) = 



Aq{x), 
Aq{~x), 



X > 
X < 0. 



(1) 

(2) 



E and E are accessible through hard exclusive electro- 
production reactions only and are new leading-twist 
functions. Similarly, in DIS, which corresponds to the 
limit ^ 0, the region — ^ < a: < ^ is absent. In this 
region the GPDs behave like meson distribution ampli- 
tudes and contain completely new information about nu- 
cleon structure. At finite momentum transfer, the first 
moments of the GPDs are related to the elastic form 
factors of the nucleon through model-independent sum 
rules. By integrating over x one can obtain for a partic- 



where 



dvcsTbh 



DVCs'^BH 



(9) 



is the interference term. Tbh is real to the lowest order 
in the QED fine structure constant. The lepton energy 
fraction y and Bjorken variable xb are defined as: 



y = 



p_q 
p ■ k^ 



Xb = 



2p-q 



(10) 



In the notation of Ref. 3 e = 2xbM/Q, where M is the 
nucleon mass, (j) is the angle between the leptonic plane 
and the proton-photon production plane. 

Contraction of the leptonic and hadronic tensors gen- 
erates an azimuthal angular dependence of each of the 
three terms in Eq.® In a frame with the z-axis 
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DVCS 



Bethe-Heitler 



FIG. 1: (Color) Feynman diagrams for DVCS and 
Bethe-Heitler processes contributing to the amplitude of 
ep ep7 scattering. 



along the virtual photon, the dependence of the ampli- 



J 



tudes on (f) yields a finite sum of Fourier harmonics. The 
amplitude of deeply virtual production of a photon has 
been derived up to twist-3 accuracy Q. In the nota- 
tion of Ref. Q , the helicity-dependent angular moments 
are presented in a series of sin(n0), with n — 1,2,3. 
Only n = 1 is a twist-2 quark matrix element. The 
n — 2 terms are twist-3 and the n — 3 terms are twist-2 
double-helicity-flip gluon transversity terms, which are 
kinematically suppressed. If these terms are omitted, 
then at the twist-2 level, the BH, DVCS, and interfer- 
ence contributions to the total cross section in Eq.® 
read 



l^vcs 

I 



1 



cos ( 



cf " cos 20} 



XBy3A2Pi(0)P2(<^) 



{cq + Ci cos + Si sir 



(11) 

(12) 
(13) 



where Vi{(j)) and 7^2 (<^) are the BH propagators. 

Experimentally, the simplest observable to measure 
GPDs is the beam spin asymmetry (Aljj)- The largest 
contribution to this observable arises from the imaginary 
part of the interference of the DVCS and the BH ampli- 



where K is the \/~t/Q^ power-suppressed kinematical 
factor and A is the beam polarization. Thus, at the 
twist-2 level, the helicity-dependent cross-section differ- 
ence will be: 



tudes, Si 



is the most interesting Fourier harmonic 



since it is linear in the Compton Form Factor (CFF) [Ic 



(14) 



(0) . (a^ 



-16Ky{2~y) 
XBy^A^Vi{cf>)V2{cl>) 



xb 



xb 



-{Fi+F2)H 



AR'P 



F2£) sin ( 



(15) 



r 



Measuring the imaginary part of the CFF \H,Ti,,£] 
gives access to the GPDs H , E] at the specific kine- 
matical point X ~ S^. From a phenomenological point 
of view, the extraction of the GPDs from data requires 
an extensive experimental program and detailed analy- 
sis with controlled theoretical corrections. This field is 
rapidly expanding (see e.g. Refs. [ll|, [H, [l3|)i and as 



a first step, the phenomenological parameterization of 
GPDs is used to fit the experimental data. 

The azimuthal dependence of the beam spin asym- 
metry has been measured by the CLAS [3] and HER- 
MES 0] Collaborations with electron (4.25 GeV) and 
positron (27.6 GeV) beams, respectively. The longitudi- 
nal target-spin asymmetry has been reported by CLAS 
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[l3| • The dependence of the hehcity-dependent and 
hehcity-independent cross sections was measured by the 
Hah A Collaboration at Jefferson Lab [l8[. The beam 
charge asymmetry [l9j and transverse target spin asym- 
metry [201 have also been reported by HERMES Collab- 
oration. Recently, CLAS Collaboration published mea- 
surements of the DVCS beam spin asymmetry over a 
wide kinematical range using high statistics data ob- 
tained with a 5.77 GeV longitudinally polarized electron 
beam [2l[. In this paper, new results on the beam spin 
asymmetry in DVCS at 4.8 GeV using the CLAS detec- 
tor are presented. Measuring the same observables at 
the same (xB^ Q^, t) but different energies, will yield 
different combinations of the terms of interest, and thus 
contribute to the determination of GPDs. In this analy- 
sis the number of single photon events has been extracted 
using the method developed in Rcf. 15]. A fit to the line 
shape of the missing mass squared distributions in the 
reaction ep — > epA was employed in each kinematical 
bin. The and t dependences of the sin moments 
have been extracted and are compared with theoretical 
predictions. 



II. EXPERIMENT 



The measurement was carried out using the CLAS de- 
tector [251 in Hall B at the Thomas Jefferson National 
Accelerator Facility. The CLAS detector system is based 
on a toroidal magnet that consists of six superconducting 
coils. The coils are positioned symmetrically around the 
beam line and provide a toroidal field in the azimuthal 
direction. Each of the gaps between the coils is instru- 
mented with an identical detector package covering typ- 
ically 80% of 27r in the azimuthal direction. 

The CLAS detector package includes three regions of 
Drift Chambers (DC) and a set of Scintillator Coun- 
ters (SC) for the tracking and identification of long-lived 
charged particles using the time-of-fiight technique. The 
DC and SC cover the laboratory polar angular range 
from 8° to 140°. In the forward region, CLAS has gas 
threshold Cerenkov Counters (CC) and Electromagnetic 
Calorimeters (EC) for electron identification. The elec- 
tromagnetic calorimeters are also used for detection and 
identification of photons and neutrons. The momentum 
resolution of the system is approximately 1%. The effi- 
ciency of single charged particle detection is > 95% in 
the fiducial volume of the detector. 

The data used in this analysis were taken during 
February-March of 2000. The 4.8 GeV longitudinally 
polarized electron beam was incident on a 5-cm long 
liquid- hydrogen target. The average beam current was 5 



uA. The total integrated luminosity was 1.28 fb^^. The 
CLAS Data Acquisition system (DAQ) was triggered by 
a coincidence of signals from the Cerenkov Counters and 
the Electromagnetic Calorimeters. The DAQ rate was 
1.1 kHz at 95% live time. The beam polarization was 
measured several times during the run using the Hall B 
MoUer polarimeter. The average polarization was 70% 
with an uncertainty of ±3%. 

A total of 1.26 M triggers have been processed, of 
which about 18% contained a properly identified elec- 
tron. For the physics analysis, events with only one de- 
tected electron and one detected proton were used. For 
systematic checks, the final states (ep7) and (ep77) were 
analyzed. 

III. PARTICLE IDENTIFICATION 



At the initial stage of data processing, the three- 
momenta of all final state particles are defined and 
preliminary particle identification (FID) is performed. 
Later, in the physics analysis, these PIDs are refined us- 
ing the knowledge of the event kinematics. 

For the identification of electrons, the energy de- 
posited in the EC and the number of photoelectrons de- 
tected in the CC are used. This method relies on the cor- 
rect reconstruction of the shower energy and on the anal- 
ysis of the electromagnetic shower profile. Fiducial cuts 
in the calorimeter plane are used to reject particles that 
pass close to the edges and lose part of their energy out- 
side of the calorimeter detection volume. In Fig. [5] the 
distribution of the number of photoelectrons detected in 
the CC, iVp.e., vs. the ratio AEec/P is shown after the 
fiducial cuts have been applied. Here AEec is the elec- 
tromagnetic shower energy detected in the EC and P is 
the momentum of the particle reconstructed in the DC. 
The lower horizontal band corresponds to negative pions. 
The vertical band, centered at AEec/P — 0.27 (aver- 
age electromagnetic shower energy sampling fraction of 
the EC), corresponds to electrons. A cut Np,e. > 2.5 is 
used to reject most of the pions. Then, the distribution 
AEec/P is fitted by a Gaussian function and a ±3ct cut 
around the mean is imposed to further clean the electron 
sample. 

Momentum and time-of-fiight (TOE) analysis are per- 
formed to identify protons. The CLAS TOF system al- 
lows the separation of protons from kaons up to 2 GeV/c 
and protons from pions up to 3 GeV/c momentum. In 
the kinematics of this analysis, protons have momenta 
smaller than 1 GeV/c. Protons are selected using a cut 
on the vertex time. At, calculated by: 

At = t- R/P~t,. (16) 
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FIG. 2: (Color) Distribution of the number of photoelectrons 
detected in the CC vs. the sampling fraction of the EC. 



Here, t is the time measured in the SC and R is the 
track length from the production vertex to the scintil- 
lator plane, which is determined by the tracking rou- 
tines. In Eq. (fTH]). f3 = p/y/p^ + m^, where p and m 
are the proton momentum and rest mass. The event 
start time, tg, is defined by the electron in the event as 
ts = — Re/c, where t,, is the time measured in the SC 
for the electron track and Re is the electron track length 
from the production vertex to the scintillator plane. In 
Fig. [3l the dependence of At on the momentum is pre- 
sented for positively charged particles assuming they are 
protons. The horizontal band centered at ns corre- 
sponds to protons. The lower band, emerging from neg- 
ative At's and closing to the proton band at high mo- 
mentum, corresponds to positive pions. A ±1 ns cut, 
shown by the dashed lines, is used to select protons. 

While final beam spin asymmetries for DVCS-BH 
events were obtained using fits to the missing mass dis- 
tributions in the reaction ep — > epX, events with epj 
and ep77 are used to extract the fit parameters and to 
study systematic uncertainties. The neutral particles in 
CLAS are detected and identified in the electromagnetic 
calorimeter. Neutrons and photons in the EC are sepa- 
rated by analyzing the speed of the neutral hits (Pec)- 
For the selection of photons, a cut on (3ec > 0.85 is used, 
see Fig. The energy of the photon is reconstructed 




P (GeV/c) 

FIG. 3: (Color) Vertex time distribution as a function of mo- 
mentum for positive tracks. The particle velocity is deduced 
from the momentum using the proton mass. 

using the energy deposited in the EC, corrected for the 
sampling fraction of the calorimeter of 0.27. 



IV. KINEMATIC CORRECTIONS 



Two types of corrections have been applied to the mea- 
sured momenta of electrons and protons. The first one 
is a correction to the proton momentum to account for 
energy loss in the material of the target and in CLAS. 
The second one is a correction to account for small un- 
certainties in the magnetic field map and in the drift 
chamber alignment. 

The energy loss corrections for protons are derived us- 
ing the GEANT simulation package for the CLAS de- 
tector. Protons are simulated in the kinematic region 
covered by the experimental data. The ratio of the gen- 
erated over the reconstructed momenta is parameterized 
as a function of the reconstructed momentum using a 
polynomial function. This function is then used to cor- 
rect the measured momentum in the analysis of the ex- 
perimental data. 
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FIG. 4: (Color) Distribution of the energy measured in the 
EC vs. /3 for neutral particles detected in the EC. The vertical 
band at /3 = 1 corresponds to photons. 



To correct for the effects of the DC misalignments and 
the uncertainties in the magnetic field distribution, two 
methods, applied to kinematically complete reactions, 
were studied. The first method corrects only the magni- 
tude of the electron and proton momenta. In the second 
method, a complex fitting algorithm was used to derive 
corrections for momenta and angles of all charged par- 
ticles in an exclusive event. These corrections depend 
on the momentum, angle, and charge of the final state 
particle and do not depend on the particle type. 

In the first method, the electron momentum correc- 
tions are derived using the Bethe-Heitler (BH) events 
associated with radiation of an energetic photon by the 
incoming electron. The kinematics of the secondary elec- 
trons in such events are similar to the kinematics in 
the deep inelastic scattering regime. For the selection 
of the BH events, cuts on the missing mass squared, 
on the transverse component of the missing momentum, 
and on the difference of azimuthal angles of the electron 
{(j)e) and the proton (t/fp) in the reaction ep —^ epX are 
used. In Fig. [5] a) the missing transverse momentum 
distribution is shown for events with the missing mass 
squared within ±0.1 (GeV/c^)^. The corresponding dis- 
tribution for the azimuthal angle difference, \(f>e — (t>p\, 
is shown in Fig. [5] b) with a dashed line histogram. A 



cut V(-Pm)^ + (-Pm)V-F'm < 0.01 has bccu applied to se- 
lect events with missing momentum in the direction of 
the beam (the radiated photon is in the direction of the 
incoming electron). Here Pm is the magnitude and P,f, 
and are the transverse components of the missing 
momentum. In Fig. [5] b) , the solid line histogram cor- 
responds to the azimuthal angle difference after the cut 
on the transverse component of the missing momentum. 
The distribution peaks at 180°, as it should for the elas- 
tic scattering events. A cut 178° < |0e — <t>p\ < 182° was 
used to select the final event sample. 

The scattered electron momentum, Pj, was calculated 
using the measured polar angles: 

Pf 



Pf = 



Pf 



1 + 

M 



1 



cos ( 



— cos 9e) 
(cos 9e + COS Or 



sin 9p 



sm( 



1), (17) 



where 9e and 9p are the polar angles of the electron and 
the proton, respectively, and M is the proton mass. Pf 
is the energy of the interacting electron after radiation. 
The ratio of over the measured momentum was pa- 
rameterized as a function of the measured angles and 
momenta, and used as a correction factor to the mea- 
sured electron momentum. 

For the proton momentum, corrections derived for the 
TT"*" from the reaction ep emr~^ with a missing neutron 
are used. The en^n events are selected using a ±3ct cut 
around the neutron mass peak. The pion momentum 
was calculated after applying corrections to the electron 
momentum and from the assumption that the missing 
particle is a neutron. As above, the ratio of the cal- 
culated and measured momenta was parameterized as a 
function of the measured angles and momenta, and was 
used as a correction factor. 

In the second method, we used the Bethe-Heitler and 
7r° production events. The BH events are selected in the 
same way as above. In order to identify the pion events, 
two photons, in addition to an electron and a proton, 
have to be detected. A cut on the invariant mass of the 
two photons was used for the final selection of 7r° events. 
Using the missing mass squared of the photon and tt" as 
a constraint in the fit, corrections for the momenta and 
polar angles of the electron and the proton are derived. 

Both methods give comparable results. For the final 
analysis, corrections obtained by the second method are 
used. In Fig. [Sla) the missing mass squared distribution 
for BH events is shown before (dashed line histogram) 
and after (solid line histogram) corrections are applied. 
After corrections are applied the position of the miss- 
ing mass peak has moved to the correct value and the 
width of the peak has improved. In Fig. [6] b) the im- 
provement on the mean values of the missing mass (Mj) 
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FIG. 5: (Color) Selection of the BH events, a) The distribu- 
tion of the X and y components of the missing momentum. 
The strong peak at zero corresponds to BH events, b) The 
electron and proton azimuthal angle difference before (dashed 
line histogram) and after (solid line histogram) the cut on the 
X and y components of the missing momentum. 



distributions of BH events as a function of CLAS sector 
are shown. The variation of the centroids is less than 
0.001 GeV^ after momentum corrections. 
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FIG. 6: (Color) a) The missing mass squared of ep epX for 
the BH events before (the dashed line histogram) and after 
(the solid line histogram) momentum corrections are applied, 
b) The dependence of the Gaussian means of the missing mass 
squared distributions on the CLAS sector before (asterisks) 
and after (crosses) the momentum corrections. 
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DVCS ANALYSIS 



The beam spin asymmetry in DVCS is studied in the 
kinematical region of W > 2 GeV and 1.4 (GeV/c)^ < 
< 2.8 (GeV/c)^ using the reaction ep epXwheie 
only one electron and one proton are detected in CLAS. 
The kinematical coverage is shown in Fig. [71 

The azimuthal angular modulations of the beam spin 
asymmetry (Am) were extracted for three bins in Q^, 
shown with horizontal lines in the left graph of Fig. [3 
for the transferred momentum range from 0.1 (GcV/c)^ 
to 0.4 (GeV/c)^. The t dependence of the angular mod- 
ulations have been studied in three bins of t, shown by 
the vertical lines in the right graph, for the range 
from 1.4 (GeV/c)2 to 2.5 (GeV/c)^. 

For each and t bin, the data were divided into 11 
bins of (/>p^, the azimuthal angle of the proton-photon 
production plane to the electron scattering plane (see 
Fig. [5]). Note that in our notation is defined as 
(j) — TT ~ (j)p~f. In each kinematical bin, the number of 
ep ep7 events for positive (A''+) and negative (A'^^) 
helicity states of the beam electrons and for the helic- 
ity sum {N^) have been extracted using a fit to the line 
shape of the missing mass squared distributions. The 
beam spin asymmetry was calculated as: 



1 N+ - 

A ^ 7 7 

^LU = — 



P. 



(18) 



where Pg is the average polarization of the beam. The 
helicity-related charge asymmetry has been measured us- 
ing the inclusive electron yield. It was found to be 0.7% 
and is included in the estimation of the systematic un- 
certainties. 

After omitting terms suppressed by an order of —t/Q^ 
or higher in Eqs.lTl]), and Alu can be 

written as: 



A 



a sm ( 



LU 



1 + J cos ( 
where the coefficients a' and 7 are: 



7 



y 



(19) 

(20) 
(21) 



In our analysis, the sin(j)p-y moment of Am, which 
is a linear combination of GPDs (see Eg fTS]) . has been 
extracted in each and t bin using fits to the (f>p.y de- 
pendence of Am with the function: 



In this representation all the higher-order terms in the 
azimuthal angle have been combined in the parameter (3. 
Note that power-sugpressed contributions could modify 
Eqs. dUl) and dUl) 



A. The Missing Mass Technique 



The main challenge in the DVCS analysis using the re- 
action ep — > epAis the separation of the single photon 
events from the more than one photon events, mostly 
tt'^. In the kinematics of deep inelastic scattering, the 
CLAS resolution of the missing mass squared is not suf- 
ficient to separate these two final states cleanly. For 
the separation of epj and 6^77 events, a fit to the line 
shape of the missing mass squared distribution is used. 
This technique for the extraction of the DVCS infor- 
mation in the reaction ep — + epX has already been 
used in the DVCS analysis of the CLAS data at 4.2 GeV 
beam energy J^] . The high statistics of this data set and 
the comprehensive simulation, enable detailed studies of 
the systematic uncertainties associated with the missing 
mass technique. 

The main contributions to the missing mass squared 
distribution in the range -0.1 GeV/c^< < 0.2 

GeV/c^ come from three processes: 1) single photon 
production; 2) 7r° production; and 3) the radiative pro- 
cesses associated with photon and 7r° production, such 
as the ep^R-f or ep^Rir^ final states. Here 77? is the 
photon radiated by the incoming or outgoing electron. 
There are also background events associated with parti- 
cle misidentification and multi-pion production. In the 
fit, it is assumed that the missing mass squared distribu- 
tions corresponding to single photon and single tt*^ pro- 
duction have a Gaussian shape defined by the detector 
resolution. The radiative tails associated with 7 and tt'^ 
production, and the background from particle misidenti- 
fication and multi-pion production, is fitted by a polyno- 
mial function. A sum of two Gaussians and a polynomial 
function is used in the final fit. 



F = ■ + N^o ■ G^o + P3 • P0I3 



(23) 



A 



LU 



a sin if)pj + P sin 2(/)„ 



(22) 



where the fit parameters A'^ and N^o are the number 
of single photon and single pion events, respectively. 
The fit parameter P3 is the relative magnitude of the 
background. The shape of the background, Pol^, has 
been determined by a fit to the tails of the missing mass 
squared distribution. The mean values and the standard 
deviations of the Gaussian functions G-y and G^^o are de- 
termined by the fits to the missing mass squared dis- 
tributions of the BH, ep — > e'p^R, and the single pion, 
ep e'p7r°, events. The identification of the BH and 
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FIG. 7: (Color) The kinematical ( 



age of the data set. (Left) The vs. xb distribution is shown with our bin divisions 



and (right) the vs. t distribution is shown. The boxes indicate the kinematical bins used for this analysis (see Tables Hill 
and TTh. 




FIG. 8: (Color) The kinematics of electroproduction in the 
target rest frame. The azimuthal angle (pp-y is the angle be- 
tween the proton-photon production plane and the electron 
scattering plane. The azimuthal angle defined in Ref. 3 and 
presented in Eq.((T5} is = tt — <pp^. 



7r° events is conducted the same way as was described 



parameters and (Tx are the mean value and the stan- 
dard deviation of the Gaussian function. The function 
P0I3, as in Ea. ((23)) . is used to describe the radiative tail 
and the background under the missing mass peak. An 
example of a fit to the missing mass squared distribu- 
tion of ep epX for the BH events in the kinemati- 
cal region 1.0 (GeV/c)2< < 1.35 {GeY/cf and 0.1 
(GeV/c)2< -t < 0.4 (GeV/c)2 is shown in FigHi. A fit 
to the distribution of the single pion events for the 
same kinematical bin is shown in Fig. [5Jd. The parame- 
ters AIx and cTx have been determined for each and t 
bin. 

The missing mass squared distributions for each he- 
licity state and for the helicity sum in each kinematical 
bin are fitted using the function in Ea. (|23p . An exam- 
ple of a fit for a typical bin, (j)pj from 35° to 70°, 1.0 
(GeV/c)2< q2 ^ ;l 35 (GeV/c)^ and 0.1 (GeV/c)2< 
-t < 0.4 (GeV/c)2 is shown in Fig. [TOl First, a fit with 
a third order polynomial is performed to the points out- 
side of the peak region, -0.1 (GeV/c)2< < -0.08 



above. For the extraction of the Gaussian parameters, (GeV/c) and 0.12 (GeV/c) < < 0.2 (GeV/c) , to 



the function f = A ■ G. 



7(Tr°) 



P0I3 is used to 



fit the missing mass squared distributions of identified 
BH and tt^ final states individually. Here, the fit param- 
eter A is the amplitude of the Gaussian function. The 



determine the shape of the polynomial P0I3, see Fig.lTUk. 

Then, using the corresponding mean values and stan- 
dard deviations of the Gaussian functions Gj and G^o, 
fits to the A/J distributions for the positive (Fig. [TUb) 
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FIG. 9: (Color) Fits to identified epj a) and epiT° b) final states to determine the mean values and standard deviations of 
the Gaussian functions G-y and G,r°- The solid line is the fit function, the dashed line is the fitted Gaussian function, and the 
dashed-dotted line is the polynomial fit to the radiative tail and the background. 



and negative fFig. \Wb) beam helicity states, and for the 
helicity sum (Fig. [TUti) were performed. In Fig. [TUb . 
Fig.fTOb. and Fig.fTOH the solid lines are the fit functions 
as defined in Ea. (P5)) . the dashed hues are the Gaussian 
functions for the single photon events, the dashed-dotted 
lines correspond to the Gaussian functions for the single 
7r° events, and the dashed-dot-dot-dot lines correspond 
to the polynomial background. 

Fits to the distributions have been performed in 
11 (ppj bins for each and t bin. The extracted num- 
ber of photon events are used to calculate Aljj for each 
kinematical point. 



gether and the missing mass squared distributions of the 
mixed samples are fitted using the function in Eq.(P5|) 
to reconstruct the number of photon and pion events. 
The difference between the initial and the reconstructed 
number of photon events is taken as a measure of the un- 
certainty. Using the fit results to the distributions 
for the DVCS analysis as a guide, the mixed samples with 
different statistics and the relative ratios of the photon 
and pion events have been studied. 



A. Studies With Experimental Data 



VI. STUDIES OF THE ANALYSIS METHOD 



The mean values of the missing mass squared distri- 
butions for the single photon and 7r° final states are sep- 
arated by approximately one-a of the fitted Gaussian 
distributions. This implies that any small systematic 
uncertainties in the energy calibration will directly af- 
fect the fit results. The accuracy of the determination 
of iV^ has been studied using real and simulated data. 
In both cases, single photon and tt*^ events are mixed to- 



The main advantages of using measured data are the 
correct representation of the background shape and the 
missing mass squared resolutions. First, separate sets 
of identified epj and epn^ final states have been created. 
The method of identification of the photon and pion final 
states has been explained before. The initial number 
of events in each final state for the mixing in a given 
kinematical bin is determined by using a single Gaussian 
and a polynomial fit to the distributions. Then the 
two samples were mixed and the number of events in each 
final state in the mixed sample was then reconstructed 
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FIG. 10: (Color) Example of fits obtained using Eq. ([23} for one bin. a) Fit to the background using the end points of 
the AfJ distribution, b), c), and d) are the fits to the distributions corresponding to positive, negative, and summed beam 
helicity states, respectively. In b), c), and d) the fit parameters are the number of single photon and n'^ events. 



by a fit to the line shape of the missing mass squared 
distribution by using the function in Eg. ((23)) . 

Two cases of event mixing are considered. In the first 
case, the mixing ratio of the ep7 and epir^ events is 
kept constant at 1.5, while the total number of events 
in the mixed distributions was changed. The red points 
in Fig. [TTk show the ratio of the initial number of pho- 
ton events to the number of reconstructed events (from 
a fit to the distribution of mixed events). The num- 
bers under each point on the plot correspond to the ini- 
tial number of epj events. Within statistical uncertain- 
ties, the reconstructed and the initial number of photon 



events are the same, although there is a general trend 
of the reconstructed numbers to be a few percent larger 
than the initial number. 

In the second case, the number of ep7r° events was 
kept constant at 360, while the number of epj events in- 
creased with each trial. In Fig. [Tib the ratio of the initial 
to the reconstructed number of single photon events is 
shown with red points. The numbers under each point 
show the ratio of the epj to epn^ events in the mixed 
sample. Again, within statistical uncertainties, the ra- 
tios of the initial to reconstructed number of single pho- 
ton events are consistent with 1. 
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FIG. 11: (Color) Ratios of initial to reconstructed single photon events a) for a constant N^/N^o ratio and b) for a constant 
number of events in the mixed distribution. In a) the numbers under each point on the plot correspond to the initial number 
of ep7 events. In b) the numbers under each point show the ratios of the epj to epTr" events in the mixed sample. 



The studies with mixed samples of identified photon 
and pion final states using the experimental data show 
that the number of single photon events can be recon- 
structed using the fit to the distribution with better 
than 5% accuracy. 



B. Studies With Simulated Data 



Another set of tests has been performed using a realis- 
tic event generator psj for the ep — > epj and ep — s- epir^ 
reactions. The response of the CLAS detector is sim- 
ulated using a GEANT-based program. The simulated 
data have been processed with the same event recon- 
struction and physics analysis algorithms that have been 
used for the measured data. Events from both samples 
(that passed all analysis cuts) in a given kinematical bin 
are mixed together and fits to the line shape of the 
distributions of the mixed events have been performed 
to extract the number of single photon and single 7r° fi- 
nal states. The parameters for Gj and Gt^o are derived 
from the fits to the missing mass squared distributions 
of the ep epj and ep —> epTr" events before mixing, 
respectively. 

Examples of the fits to the Mj distributions for the 
different helicity states of the simulated data are shown 
in Fig. 1121 The open circles with error bars correspond 



to the mixed event distributions. The solid line is the 
fit using the function in Eq. (P5)l to these points. The 
open squares represent the distributions of the initial 
photon events. The dashed lines are the resulting G^ 
from the fit to the mixed distributions. The asterisks 
and the dashed-dotted lines are the same for the single 
pion final state. The dashed-dot-dot-dotted line (shown 
only in Fig. [T^b) corresponds to the fitted polynomial 
background. 

In Fig. [13] the ratio of the number of single photon 
events, obtained from fits to the missing mass squared 
distributions of the mixed event sample, to the number of 
reconstructed photon events that passed all analysis cuts 
are shown for all ipp^ bins at 1.0 (GeV/c)^< < 1.35 
(GeV/c)2 and 0.1 (GeV/c)2< t < 0.4 (GeV/c)^. The 
crosses and the circles represent the ratios for different 
helicities. The asterisks show the ratio of these ratios. 
The extracted number of events by the fit is always lower 
than number of events in the mixed sample by about 
20%. This is due to the fact that there is no background 
and radiative effects in the simulated data, while the 
fit procedure finds a non-zero contribution for the poly- 
nomial function at the expense of the photon and pion 
events. The important information is indicated by the 
asterisks, which show that for each ^p-^ bin, these ratios 
are the same for both helicities. 

In Fig. I14[ the beam spin asymmetry, as defined in 
Eq. fTB)) . is presented as a function of (f>p^. The simulated 
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FIG. 12: (Color) Simulated missing mass squared distributions for positive a) and negative b) helicity for one (/f>p-y, Q^, and t 
bin. The curves are results of fits using the function in Ea. (|23p . In the simulation radiative effects are not included. 



(f)pj dependence of the asymmetry is shown with circles. 
The crosses represent the asymmetry obtained using the 
fit procedure. The soHd Hne is a fit using the function in 
Eg. ((22)) to the extracted asymmetry. The dotted line is 
a fit to the simulated asymmetry. The sin(/)p-y moment 
for the extracted asymmetry is approximately 3% higher 
than the simulated moment. This difference is used as 
the systematic uncertainty in the determination of the 
moments of the azimuthal asymmetry in DVCS. 

C. Dependences on the Parameters of G-j and G^o 



within a ±1 RMS range. Thus, a total of 5^ = 625 dif- 
ferent sets of means and standard deviations were tested 
for each t and bin. 

In Fig. 116b . the distribution of the sin (pp^ moments ob- 
tained for different sets of Gaussian parameters is shown 
for the Q2 and t bin of 1.35 to 1.7 (GeV/c)^ and 0.1 
to 0.4 (GeV/c)2, respectively. The RMS of this distri- 
bution is 0.015. The same distribution for the sin20j,^ 
moment has an RMS of 0.013, see Fig. Mb- These RMS 
values were taken as the systematic uncertainties of the 
moments due to the kinematical dependence of the pa- 
rameters of G-y and G^o within the kinematic bin. 



Additional sources of systematic uncertainties are the 
dependences of the mean values and standard devia- 
tions of the functions G-y and G^^o on the kinematics 
of an event. In Fig. [151 distributions of these parame- 
ters obtained in smaller sub-bins of t and within the 
g2 bin of 1.7 (GeV/c)2 to 2.8 (GeV/c)^ and t bin of 
0.1 (GeV/c)2 to 0.8 (GeV/c)^ are shown. Single pho- 
ton [Bethe-Heitler events, selected as shown in Fig. [S] 
events in this t and bin were divided into 30 sub-bins 
fFig. ll5l (a.b)). while single pion events were divided into 
21 sub-bins (Fig. [TSTcd)). Physics asymmetries in each 
t and bin were extracted using the central values of 
these parameters in the given bin. To estimate the sys- 
tematic uncertainties, the whole fitting procedure was 
repeated with additional four values of these parameters 
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FIG. 13: (Color) The (j>pj dependence of ratios of the number 
of single photon events extracted using the fits to the missing 
mass squared distribution of the mixed events to the number 
of single photon events that passed the analysis cuts for the 
positive (open circles) and for the negative (crosses) beam 
helicity states. The asterisks represent the ratio of the two 
ratios. 



FIG. 14: (Color) The simulated and extracted Alu- The 
crosses are the extracted asymmetries and the solid line is the 
fit to these points using the function in Ea. (|22p . The open 
circles and the dotted line are the same for the simulated 
asymmetry. 



16 




-0.01 0.01 0.02 
Mean (GeV^) 




0.02 0.04 

X 



Sigma (GeV^) 




0.02 0.04 
Mean M^ (GeV^) 




0.02 0.04 
Sigma M^ (GeV^) 



FIG. 15: Distribution of the means and standard deviations of the Gaussian functions for the photon and pion peaks for 
sub-bins of one kinematical bin: 1.7 [GeV/cf < < 2.8 (GeV/c)^ and 0.1 (GeV/cf < \t\ < 0.8 (GeV/cf 
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FIG. 16: Distribution of the sin (cl) and sin '2(f)p-^ (b) 
moments of the beam spin asymmetries for the bin of 
1.35 (GeV/c)^ to 1.7 (GeV/c)^ and 0.1 (GeV/c)^ < |t| < 
0.4 (GeV/c)^, extracted from the fits to the missing mass 
squared distributions using different values of means and 

within ±1 RMS. 



18 



VII. RESULTS 



Our results on the beam spin asymmetry, calculated 
using Eq. , are presented in Tables U and |TT1 In or- 
der to extract sin^p.^ (a) moments as a function of 
and the azimuthal angular dependence of Am was fit- 
ted with the function presented in Eq.([22|l. In Fig. [17] 
the (jip-y dependences of Am for three bins of Q'^ for 
0.1 (GeV/c)2 <-t< 0.4 (GeV/c)2 are shown. The az- 
imuthal angular dependence oi A^jj for three bins of t at 
1.4 (GeV/c)2 < Q2 < 2.5 {GeY/cf is shown in Fig.[THl 
The solid line is a fit of the extracted asymmetry using 
the function presented in Eq. (P^ . In both figures only 
statistical uncertainties are shown. 



< —t > 


0.19 [{GeV/cf 


0.30 [{GeV/cf 


0.46 [(GeV/c)^] 


< (/) > 


At IT 


^ A T TT 


A T TT 


^ At n 


At TT 


^ At TT 


1 5° 


0.0014 


0.0102 


-D 0009 


0.0071 


0.0141 


0.0049 


6.5° 


0.0449 


0.0219 


0.0024 


0.0185 


0.0209 


0.0131 


22.5° 


0.0298 


0.0260 


0.1143 


0.0266 


0.0969 


0.0204 


52.5° 


0.1354 


0.0402 


0.1871 


0.0503 


0.1998 


0.0450 


102.5° 


0.2361 


0.0503 


0.1191 


0.0906 


0.2461 


0.0984 


180.0° 


-0.0858 


0.1092 


0.1608 


0.1994 


0.0629 


0.1908 


257.5° 


-0.2164 


0.0533 


-0.0513 


0.1210 


-0.3280 


0.1797 


307.5° 


-0.2105 


0.0389 


-0.1756 


0.0548 


-0.1679 


0.0515 


337.5° 


-0.0691 


0.0253 


-0.0594 


0.0260 


-0.0653 


0.0194 


353.5° 


0.0111 


0.0218 


0.0010 


0.0179 


-0.0377 


0.0129 


358.5° 


0.0075 


0.0110 


-0.0035 


0.0076 


0.0045 


0.0053 



<Q''> 


1.22 [{GeV/cf 


1.51 [{GeV/cf 


2.04 [{GeV/cf] 


< (t>> 


Alu 


SAlu 


Alu 


SAlu 


Alu 


SAlu 


1.5° 


-0.0051 


0.0056 


0.0007 


0.0055 


0.0021 


0.0060 


6.5° 


0.0289 


0.0139 


-0.0005 


0.0134 


0.0315 


0.0154 


22.5° 


0.0794 


0.0221 


0.0660 


0.0196 


0.1054 


0.0217 


52.5° 


0.1427 


0.0381 


0.1826 


0.0322 


0.1544 


0.0454 


102.5° 


0.2028 


0.0589 


0.2080 


0.0475 


0.2237 


0.0750 


180.0° 


0.0227 


0.1062 


0.0703 


0.1146 


0.0071 


0.1307 


257.5° 


-0.1394 


0.0582 


-0.2592 


0.0539 


-0.1471 


0.0906 


307.5° 


-0.1662 


0.0431 


-0.2173 


0.0337 


-0.1451 


0.0463 


337.5° 


-0.0862 


0.0203 


-0.0631 


0.0185 


-0.0643 


0.0220 


353.5° 


-0.0147 


0.0135 


0.0130 


0.0129 


-0.0308 


0.0158 


358.5° 


-0.0032 


0.0063 


-0.0001 


0.0058 


0.0006 


0.0063 



TABLE II: Calculated beam spin asymmetry, Alu, at each 
(jjp-y for three t bins, see Fig. [T] The average values of a;s and 
are shown in Table HVl 



TABLE I: Calculated beam spin asymmetry, Alu, at each ^ 
(f)pj for three bins, see Fig. [7] The average values oi xb 
and t axe shown in Table Hill 

The systematic uncertainties on a have been evalu- 
ated using the studies presented in the previous section. 
Besides the systematic uncertainties due to the fits, 
3%, and the determination of the parameters of Gy and 
Gj^o, ±8%, there is an uncertainty in the calculation of 
the Alu from the determination of the beam polariza- 
tion, Pe, and an uncertainty due to the charge asymme- 
try. The uncertainty on the beam polarization was esti- 
mated to be ±3%. The helicity-dependent charge asym- 
metry was about 0.7%. The estimated total systematic 
uncertainties for the extracted moments are —8.7% and 
+8.2%. 

The final results on the moments a and /3 are pre- 
sented in Table imi and Table HVl For consistency checks, 
the (/>p^ dependence of the Alu were fitted with the func- 
tion presented in Eq. (fT9| . The obtained results on a' and 
7 are also presented in Table IIIII and Table IIVI In the 
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FIG. 17: Asymmetry (Alu) as a function of azimuthal angle 
fitted according to Ea. (|22p to extract the sin(^ and sm2(f> 
moments in three bins, see Table Hill 



tables only the statistical uncertainties are quoted. 

In Fig. [191 the and ^-dependences of the sin (j)pj 
moment (a) are shown. In the graph for the depen- 
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< Xb > 


< > [{GeV/cf] 


a 




a 


7 


1.22 


0.17 


0.23 


0.181 ±0.032 


0.099 ±0.023 


0.181 ±0.032 


-0.098 ± 0.228 


1.51 


0.20 


0.26 


0.245 ± 0.028 


-0.040 ±0.021 


0.234 ± 0.024 


0.319 ±0.195 


2.04 


0.28 


0.38 


0.192 ±0.044 


0.010 ±0.030 


0.191 ±0.045 


-0.107 ±0.288 



TABLE III: Results from the fits to the (j> dependences of Alu with the functions presented in Eqs. (I19p and (|22|l . Only 
statistical uncertainties are presented. 



<-t> [{GeY/cf] 


< Xb > 


< Q2 > [(GeV/c)2] 


a 




a' 


7 


0.19 


0.23 


1.69 


0.228 ± 0.029 


-0.034 ±0.024 


0.222 ± 0.026 


0.318 ± 0.248 


0.30 


0.23 


1.74 


0.147 ±0.055 


0.033 ±0.037 


0.161 ±0.045 


-0.223 ± 0.274 


0.46 


0.25 


1.83 


0.260 ± 0.062 


-0.022 ±0.039 


0.260 ± 0.057 


0.200 ± 0.331 



TABLE IV: Results from the fits to the (p dependences of Alu with the functions presented in Eqs. (|19|) and (|22|l . Only 
statistical uncertainties are presented. 




(degree) 

FIG. 18: Asymmetry as a function of azimuthal angle fitted 
according to Ea. (|22ll to extract the smcp and sin 20 moments 
in three t bins, see Table IIVI 



dence the previously published CLAS result on a [T^ is 
shown by an open cross. The solid and dashed-dotted 
lines are the and t-dependences of the sin (ppj mo- 



ments of the beam spin asymmetry calculated using a 
Regge trajectory exchange model for the interaction of 
the photon with the proton |,24„ 25j using two treatments 
of the unitarity cuts. The dashed-dotted curves retain 
the contribution of the poles [ISl and the contribution of 
the p-p elastic cut [25|. They take into account the cou- 
pling between the 7-p and the p-p channels and are well 
under control since they rely on known on-shell matrix 
elements. The solid curves also take into account all of 
the inelastic diffractive cuts, under the assumption that 
diffractive dissociation saturates the photo-absorption 
cross section. This treatment leads to an excellent agree- 
ment with the unpolarized cross sections, as well as the 
helicity difference of the DVCS cross sections, recently 
measured in Hall A [l8| . 

The other two curves in the graphs represent calcula- 
tions of the siii(j)pj moments with a GPD-based model 
[26j . using a Regge ansatz for the <-distribution of the 
GPD [231 . The parameterization of GPDs includes the 
D-term (to ensure the polynomiality of Mellin moments 
of GPDs [11]). Helicity-dependent cross sections were 
calculated at the twist-three level and include target 
mass corrections [2^. Calculations are done for two 
different skewedness parameter values for sea quarks, 
bsea = 1 - dashed line, and bsea — 2 - dotted line. In the 
calculations the profile parameter for the valence quark 
was bvai = 1 and the Regge slope parameter a' = 1.05 
[27} . The data favor the GPD model with the sea quark 
skewedness parameter b^ea — 2. The calculations include 
the DVCS contributions to both the numerator fEq. p^ ) 
and denominator (Eq. p^ ) of A^jj- 

In the method used for separation of the single pho- 
ton events, the radiative tails were folded into the "back- 
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FIG. 19: (Color) The (a) and t (b) dependences of the sin (pp-f moments of the DVCS asymmetry. The "cross" symbol is 
the asymmetry measured by CLAS with the 4.2 GeV data. The curves are model calculations: the solid and dashed-dotted 
lines are the calculations based on a Regge model for the photon-proton interaction [2^. The dashed and dotted lines are 
from the GPD-based model [2^. Isol]. Each curve represents different parameter settings (see text for details). In the plots the 
error bars are the combined statistical and systematical uncertainties added in quadrature. 



ground" . The radiative effects are inseparable from the 
Born level ep — > epj process. Detailed radiative cal- 
culations by Vanderhaeghen et al. [1^ showed that the 
effect on the beam spin asymmetry is less than 5%. This 
is much smaller than the uncertainties of our measure- 
ments. 

The extracted sin 2(f>pj moment or cos (f>p-y denomina- 
tor term of Am is within uncertainties consistent with 
zero. The precision of our results is not good enough 
to determine small contributions from higher twists and 
the power suppressed terms. 



VIII. SUMMARY 



The beam spin asymmetry Aj^jj in the deeply virtual 
production of real photons has been measured using a 
4.8 GeV longitudinally polarized electron beam and the 
CLAS detector. The experimental A^jj for each kine- 
matical bin was calculated using the number of single 



photon events extracted from the fits to the line shape 
of the missing mass squared (Mj) distributions. Studies 
using Bethe-Heitler and tt*^ production events have been 
performed to test the validity of the fit method. It was 
found that the systematic uncertainties on the moments 
from the fit procedure are not more than 6%. Overall 
systematic uncertainties have been estimated to be less 
than 9% of the value of a. 

The sin (j)py , sin 20p^ , and cos ipp^ moments of the az- 
imuthal angular dependence of the asymmetry have been 
extracted in three bins of and three bins of trans- 
ferred momentum t (see Tables Hill and ITV)) . The and 
t dependences of the sin (pp^ moment are compared with 
theoretical calculations using the GPD-based model for 
DVGS [2^ [s^l and the photon-proton interaction based 
on a Regge model (25| . The sensitivity of our results is 
not good enough to estimate small contributions from 
higher twists and the power suppressed terms. However, 
the data put constraints on the model parameters used to 
calculate the beam spin asymmetry. Also, the data are 
reproduced by a model based on Regge poles and unitary 
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cuts. Whether this is a consequence of the quark hadron 
duahty remains to be investigated. Clearly, more data 
in a wide range of kinematics are needed to refine the 
parameters of the models. New experiments on DVCS 
and on Deeply Virtual Meson Production will provide 
further observables in the Deeply Exclusive Production 
regime that can be used in a global fit to extract the 
GPDs. 
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